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Abstract 

The dependence of the scaling properties of the structure factor on space 
dimensionality, range of interaction, initial and final conditions, presence 
or absence of a conservation law is analysed in the framework of the large-N 
model for growth kinetics. The variety of asymptotic behaviours is quite rich, 
including standard scaling, multiscaling and a mixture of the two. The dif- 
ferent scaling properties obtained as the parameters are varied are controlled 
by a structure of fixed points with their domains of attraction. Crossovers 
arising from the competition between distinct fixed points are explicitely 
obtained. Temperature fluctuations below the critical temperature are not 
found to be irrelevant when the order parameter is conserved. The model is 
solved by integration of the equation of motion for the structure factor and 
by a renormalization group approach. 
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1 - Introduction 



In growth kinetics one deals with the relaxation to equilibrium of a system 
quenched from high to low temperature 1 . The processes of interest are those 
which exhibit scaling 2 in the asymptotic time regime. Denoting with Tj and 
Tp the initial and final temperatures, these processes can be grouped into 
three classes characterized by (Tj > T c ,Tp < T c ), (Tj > T c ,Tp = T c ) and 
(Tj = T c , Tp < T c ) where T c is the critical temperature. This subdivision 
arises from renormalization group arguments 3 whereby the temperature axis 
(Fig.l) is controlled by three fixed points at T = 0, T c , T = oo and T c is 
unstable both with respect to T = and T = oo. Such a flow diagram leads 
naturally to the three universality classes listed above whose basic processes 
are those originating and terminating in a fixed point. By far, the most 
studied among these is the phase ordering process from Tj — oo to Tp — 0. 

The reason for the continuing interest in this problem is the persistent 
lack of a full understanding of scaling which is observed both in laboratory 4 
and numerical 5 experiments. In terms of the structure factor (Fourier trans- 
form of the equal time order parameter correlation function) this asymptotic 
scaling behaviour is of the form 

C(k,t) ~ L a (t)F(kL(t)) (1.1) 

where L(t) is a characteristic length which grows in time with a power law 

L(t)~tV z . (1.2) 

A scaling pattern of this type, which we refer to as standard scaling, is com- 
pletely characterized by the pair of exponents z, a and by the scaling function 
F(x). These quantities depend to a different extent on the various elements 
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entering in the specification of the process 6 which, in addition to the classes 
(Tf,Tp) discussed above, include the space dimensionality of the system, the 
vector dimensionality of the order parameter, the presence or absence of a 
conservation law and the short or long range character of interactions. 

The purpose of this paper is to explore in detail the dependence of the 
scaling properties on the totality of these elements in the framework of the 
large-N model 7 . This, at the moment, is the only available non trivial sol- 
uble model with a structure sufficiently rich to be adequate for this kind of 
investigation. The picture which emerges in the end is quite informative and 
exposes clearly the profound difference between processes with and without 
conservation of the order parameter. 

We consider a system described by an N-component order parameter 

— # 

(j>(x) = (0i (x), 4>n(x)) and by a free energy functional of the Ginzburg- 
Landau type 



by the temperature T and by the pair of coupling constants \i = (r, g) . In 
the large-N limit (N — > oo) there is a critical temperature T c (/i) ~ — r/g 
and a phase diagram (Fig. 2) in the three dimansional parameter space (T, ji) 
with a surface of critical points separating ordered states below it from dis- 
ordered states above it. The interesting portion of this phase diagram is the 
(r < 0, g > 0) sector at or below the critical surface where scaling is to be 
expected in a quench process. As we shall see in the following, the set of 
states at T = on the g-axis plays a special role since it is located at the 
edge of both the critical surface and the ordering region below it. 




(1.3) 



where 'Hlr^] contains the long range interaction and will be specified in sec- 
tion 2. The Gibbs equilibrium states -P eg [0] ~ exp(—^7i[(p\) are parametrized 
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As anticipated above, the characterization of a process requires in the 
order: 

i) specification of the space dimensionality d 

ii) specification of the initial condition. This we do by taking an initial 
structure factor of the form 

C(M) = ^ (1.4) 

where A is a constant and the value of selects the initial state of interest: 
= corresponds to an uncorrelated initial state at infinite temperature 
(Tj = oo) while 9 = 2 corresponds to the critical point (Tj = T c ) 

iii) choice between a non conserved order parameter (NCOP) and a conserved 
order parameter (COP) 

iv) specification of the short or long range nature of the interaction 

v) specification of the final state. The interesting subsets in the equilibrium 
phase diagram are: 

[Tp — T c — 0, /ii — (r — 0, g — 0)] trivial critical state at zero temperature 
[Tp = T c > 0, fj,i — (r — 0, g — 0)] trivial critical states at finite temperature 
(T-axis) 

[Tp = T c = 0, /i 2 = (r = 0, g > 0)] non trivial critical states at zero tempera- 
ture ((/-axis) 

[Tp = T c > 0,/i 3 = (r < 0,g > 0)] non trivial critical states at finite 
temperature (critical surface) 

[Tp < T c , fx 3 = (r < 0, g > 0)] phase ordering region. 

It is convenient to regard the space dimensionality, the initial condition and 
the range of the interaction as forming, so to speak, the environment of the 
process, while the set (Tp, n) and the specification NCOP or COP as elements 
of discrimination which we will use to identify processes. 
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Solving the model analytically and by renormalization group (RG) we ar- 
rive at the following picture. The asymptotic scaling properties [z, a, F(x)] 
depend on (Tp,/x). There is a universality class, under each heading NCOP 
or COP, for each of the five regions (T F ,fx) listed above. In RG language 
this means that there are five fixed points (Tp,/x*). The flow in the param- 
eter space and therefore the extension of the universality classes depends on 
the relative stability of these fixed points. This in turn is regulated by the 
existence of critical dimensionalities which depend on the environment, i.e. 
initial condition and range of interaction. 

The deep difference between NCOP and COP emerges from how the 
scaling properties depend on the final state (T F , //). The most striking differ- 
ence is obtained for quenches inside the phase ordering region. It was found 
previously 8 that when the system is quenched to (T F = 0, // 3 ) the standard 
scaling form (1.1) holds only for NCOP, while for COP it is replaced by the 
more general multiscaling behaviour 

C(k,t) ~L a( - x \t)F(x) (1.5) 

where also the exponent a depends on x = kL(t). We find now that, 
with some modifications to be discussed below, this basic distinction NCOP- 
standard scaling and COP-multiscaling holds not just for quenches to T F = 0, 
but for quenches anywhere in the phase ordering region (T F < T c ,/i 3 ). Fur- 
thermore, while temperature perturbations with < T F < T c are irrelevant 
for NCOP, it is not so for COP. For quenches elsewhere, i.e. on the critical 
surface, standard scaling holds both for NCOP and COP. However, while 
with NCOP (Tp,/i) affects a with no impact on z, the opposite occurs for 
COP. 

The analytical tractability of the large-N model 7 ' 8 ' 9 allows to expose nicely 
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the mechanism underlying the picture outlined above and to compute in 
addition to the asymptotic properties also the crossovers induced by the 
competing fixed points. The question of the extension of the properties of 
the large-N model to finite N needs to be treated with care. We shall comment 
on this in the concluding section. 

The paper is organized as follows: in section 2 the general features of 
the large-N model are presented, section 3 is devoted to the solution of the 
model by integration of the equation of motion for the structure factor and 
in section 4 the model is analysed by RG methods. Concluding remarks are 
made in section 5. 



2 - The large-N model 

The long range part of the free energy functional (1.3) is of the form 

H LR [<p} = f d d x f d d x'4>{x) ■ V{x-x')(f{x') (2.1) 



with the large distance behaviour V(x — x') ~| x — x! \~( d+a \ Considering 
a time evolution of the order parameter governed by the time dependent 
Ginzburg-Landau model and neglecting 10,11 k 2 with respect to k u for a < 2 
the equation of motion for the Fourier transform of the order parameter in 
the large-N limit is given by 7 



d(/> a (k,t) = 
dt 



wk p+a + k p R(t) <j) a (k,t)+r] a (k,t) (2.2) 



where (a = 1,...,N), w is a coefficient originating in the small momentum 
expansion of the interaction, T is a kinetic coefficient, p = for NCOP, p = 2 
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for COP, ff(k, t) is a gaussian white noise with expectations 

<ff(k,t)> = (2.3) 

< r] a (k, t)rjp(k', t') > = 2TT F k p 5 a/3 5(k + k')8(t - t') (2.4) 

and 

R(t)=r + gS(t) (2.5) 

with S(t) =< (f>a(x,t) > which is independent of a and must be determined 
self-consistently 12 . In the following we will let a and p vary continously 
since a < 2 describes long range interactions while < p < 2 describes non 
local conservation of the order parameter 13 . 

Integrating Eq.(2.2) with a random initial condition (£(k, 0) and dropping 
the label a we find 

0(M) = 4>{k, 0)D(k,t) + f dt' V (k,t')^l (2.6) 

Jo D(k,t') 

where 

D(k, t) = exp{-T[wk p+CT t + k p Q{t)\) (2.7) 

and 

Q(t) = I dt'RCt'). (2.8) 
Jo 

From (2.6) correlation functions of arbitrary order can be obtained form- 
ing products of (f)(k, t) and averaging over both initial condition and thermal 
noise. For the average order parameter we find 

< (f)(k, t) >=< <f)(k, 0) > D(k, t) (2.9) 

— * 

which shows that if the initial state is symmetric < <f>(k, 0) >= , as we shall 

— # 

assume in the following, then < (f>(k,t) >= for all time, i.e. dynamics does 
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not brake the symmetry. The more general case of a time evolution with 
broken symmetry is outlined in Appendix I. 

The structure factor < (j)(k,t)(j)(k',t) >= C(k,t)5(k + k') is given by the 
sum of two contributions 

C(k,t) = C 1 (k,t) + C 2 (k,t) (2.10) 

where 

C l {k,t) = C{k,t = 0)D 2 {k,t) (2.11) 

C 2 (k, t) = 2Tk p T F D 2 {k, t) C dt'D~ 2 {k, t'). (2.12) 

Jo 

From (2.10-2.12) one can easily verify that the structure factor obeys the 
equation of motion 

= -2T\wk p+CT + k p R(t)}C(k, t) + ZTk p T F (2.13) 

at 

which is closed by the self-consistency condition 

/rjdh. 

Let us now analyse the final equilibrium states. Since there is no sym- 
metry breaking, it is sufficient to look at the structure factor in the limit 
t — > oo obtaining (Appendix II) from Eq. (2.13) 

C(k, oo) = , Tf , + (2n) d M 2 5(k) (2.15) 

where £ = [i?(oo)]~ 1 / CT is the equilibrium correlation length. Enforcing the 
self-consistency condition (2.14) one finds (Appendix II) that there exists a 
critical temperature 

T c = -J— (2.16) 

gB(0) 1 ; 



such that 

f t,- 1 > and M = for T F > T c , 

\ r 1 = and M 2 = M 2 (T C - T F )/T C for T F < T c [ ' 

with Mq = —r/g and 

B (0) = / t^-L oc -J— . (2.18) 

Thus, in the (Tp, /x) parameter space (Fig. 2) the critical point as r and 
g are varied ( (? must be positive in order to have a well defined theory) 
spans a surface which lies on the r < sector and separates ordered states 
underneath it, where the structure factor dispalys a Bragg peak, from the 
disordered states above it without Bragg peak. Notice that from (2.16) and 
(2.18) lim^o-Tc = 0, implying that a is the lower critical dimensionality of 
the model. 

As discussed in the Introduction, in the rest of the paper we shall be 
concerned with the solution of Eq.(2.13) with values of (Tp,/x) lying on or 
below the critical surface. 



3 - Quench processes 

From (2.10) and (1.4) the formal solution for the structure factor is given by 

C(k,t) = — e - 2 n™kp+°t+kPQ(t)] + 2 TT F k p f dt'e~ 2V{wkv+a{t ~ t,)+kPm) ~ Q{t ' ))] . 
k d Jo 

(3.1) 

We must now extract the scaling properties. 
3.1 - NCOP 
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With p = by dimensional analysis we can identify the characteristic 
length 

L(t) = (2Tt) 1/a (3.2) 

which enters (3.1) in the combination x = kL(t). It is worth to rewrite (3.1) 
with this change of variable 

C(k,t) = Ae- 2T W^re- wx ° +<jT f L°x-° f dx ' x '^ e -^-^) e -^m)-Q{t')] 
x° Jo 

(3.3) 

which shows that (3.2) is the only available choice for the caracteristic length 
and therefore necessarily z = a for any process with NCOP. In order to solve 
for Q(t) we integrate (3.3) over k and we use the definitions (2.8) and (2.14) 
to derive the equation 

1 d 

2Tdi' 

where 

r r\ d k p -ZTwk a t ML 

^=/SV=^j[ ^-e- (3.5) 

with K d = {2 d ^ 1 Ti d / 2 T(d/2))~ l , A is a monetum cutoff and A (t) is the same 
as A(t) with 6 = 0. 

T F = 

Let us now make a further restriction by considering quenches to zero 
temperature. With T F = (3.4) can be integrated obtaining 

e 2TQ(t) = e 2Trt 1 " 

'to 



_ e 2TQ(t) = re 2TQ(t) 



+ gAA(t) + 2TT F g C dt'AM - t')e 2TQ ^ (3.4) 
J o 



1 + 2TgA (J* dt'e- 2Trt 'A{t') + jT rft'e" 2 ™' A(0) 



. (3.6) 
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Next, choosing t such that AL(t ) ~ 1, from (3.5) we may approximate 

T*-*(AA(A S K d ^ L dxx d - e - l = K d [M{t)] d - e /{d-6) for t < t 
{ ) { ) ~ \ I = / °° dxx d - e - 1 e- wxa for i > t 

(3.7) 

and inserting into (3.6) finally we get 



e 2TQ(t) = Ke 2Trt 



+ gAIe 2Trt f dt'e- 2Trt 'L -\t') (3.8) 

J to 



where « = 1 + gAK d A d ~ e (l - e' 2Vrto )/{d - 6)r. 

The next step is to extract the behaviour of Q(t) from (3.8) after speci- 
fying the set of coupling constants \i = (r,g). 

fi 1 = (r = 0,g = 0) 

In this case Q(t) = . Thus from (3.3) we get 

C(k, t) = AL e F(x) (3.9) 

with 

F(x) = —. (3.10) 

Note that C(k,t) obeys the scaling form (1.1) from beginning to end over 
the entire time history of the process with a = 9 . In this case dynamics 
propagates the scaling properties of the initial condition. 

M2 = (r = 0,(7 > 0) 

Setting r = in Eq.(3.8) 

e 2TQ(t) = K+ A/a r e+a _ d _ , n 

6+a-d 1 u J 
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which gives 



with d c = + (j , L = L(t ) and given by (3.10). Now scaling is no 
more an exact property obeyed over the entire time history, but only asymp- 
totically. Furthermore, Eq.(3.12) shows that d c is a critical dimensionality, 
in the sense that different asymptotic behaviours are obtained depending on 
d > d c or d < d c . For d > d c the asymptotic behaviour is the same as for 
the quench to fii , with the denominator producing a correction to scaling 
at early times. Instead for d < d c there is a crossover time t* ~ (Ag) d - d c 
such that 

C(ht\ f L 6 F(x) for t«t* 

C ^ t )~\ L d -°F(x) fort»r. (3 - 13) 
At d = d c we find a logarithmic correction due to marginality 

cAi) ~i^k F(i) - (3 - 14) 

In order to illustrate the above results we have integrated numerically 
the equation of motion (2.13) for p — 0, a — 2, d — 3, r = 0, T F — and 
different values of g . Since the crossover involves only the exponent a it is 
convenient to discuss it in terms of S{i) . The analytical behaviour of S(t) 
is obtained integrating (3.12) over k 

jQ-d 

S ^ ~ k + aIgA[L^- d - L^]/{d c - d) ' (3 ' 15) 

The case d > d c is realized taking 9 = 0, which yields d c = 2 . The double 
logarithmic plot of the numerical solution for S(t) shows (Fig.3) that the 
power law L~ d = t~ 3 / 2 is asymptotically obeyed with deviations from it at 
early times which become smaller as g decreases, i.e. as /ii is approached 
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along the /2 2 axis. Conversely for d < d c , obtained by setting 6 = 2 which 
implies d c — 4 , the numerical solution for S(t) behaves as L^ 1 ~ t -1 / 2 for 
and as L~ 2 ~ t^ 1 for i > i* (Fig.4). 



Ms = (r < 0,(7 > 0) 



With these values of r and g the quench is made inside the region of the 
coexisting ordered phases, which is the process normally considered in the 
kinetics of phase ordering. In a process of this type usually one makes the 
distinction between the early stage of exponential growth due to the instabil- 
ity generated by r < and the late stage characterized by scaling behaviour. 
However, due to the existence of crossovers,we should expect to find a scaling 
regime also at early time with the exponents appropriate to /ii or /i 2 when 
r and g are pushed sufficiently close to ji\ or /i 2 ■ In fact for short time 
(t < l/2r|r|) Eq. (3.8) yields 

e 2r Q (t) „ e 2Trt ( K + aI 9^_ [L d c - d _ L d - d] \ (3 16) 

V d c — d J 

which implies for the structure factor 

cCk - () ~ » + ^| fr Mn/(i-/ |I|f " m (3 " 7) 

showing that in front of the exponential actually there is a prefactor identical 
to (3.12). Hence if |r| is sufficiently small one can detect a scaling regime 
identical to the one discussed for quenches to fi 2 and preceding the usual 
early time regime of exponential growth. 

Instead, when time is large we can set to zero the left hand side of (3.4) 
obtaining for the structure factor the late stage scaling behaviour 

C(k,t) ~ M$L d F(x) (3.18) 
14 



for any d . This means that for quenches inside the phase ordering region 
there is not an upper critical dimensionality such that above it one obtains 
an asymptotic behaviour with a = 9. 

The crossover structure is illustrated in Fig. 5 through the behaviour of 
S(t) computed numerically for d = 3,9 = 0,cr = 2, g = 1 and decreasing 
values of r in order to explore the influence on S(t) of the power law 
associated to the // 2 -axis. While away from /x 2 , e.g. at r = —1.0 , the 
behaviour of S(t) shows only the asymptotic scaling regime S(t) ~ Mq = 1 , 
as r is decreased and /i 2 is approached for r sufficiently small, e.g. r = —0.01, 
S(t) indeed displays at intermediate times the power law L~ d ~ t~ 3 / 2 which 
we have found previously as an asymptotic behaviour in the quenches to /i 2 . 
The case with 9 = 2 yields qualitatively similar results. 

T F >0 

For quenches to a finite final temperature < Tp < T c we restrict the 
choice of the initial condition to high temperature taking 9 = in (1.4). 

m = (r = 0,g = 0) 

For quenches to the trivial critical states at finite temperature (T-axis in 
Fig.2), setting Q{t) = in (3.3), we find 

C(k, t) = Ae~ wx ° + L a T F F (x) (3.19) 

where 

F (x) = (1 - e- wx °)/wx a . (3.20) 
Thus, for these processes the asymptotic scaling behaviour is given by L^TpF^x) 
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while the strength A of initial correlations is an irrelevant parameter with 
correction to scaling behaviour ~ L~ a . 



fi 3 = (r < 0,(7 > 0) 



The solution of Eq.(3.4) in the general case has been obtained by Newman 
and Bray 12 . For large time and for a < d < 2a they find 

e -2FQ(t) _ LU ( t j ( 3 21 ) 

with 

J e = 2a - d for T F = T c 

U ~\d for Tp < T c . (3 ' 22) 

Inserting into (3.1) we find 
C(k,t) = C(k,t ) L ^ )e _ wx „ +L°T F J Q dy(l-y)-^e— v (3.23) 

where to is some microscopic short time such that scaling holds for t > to and 
where xo = kL(t ). For Tp < T c , recalling that T c > requires d > a and 
using (3.22) the asymptotic scaling behaviour is still given by (3.18) with Tp 
acting as an irrelevant perturbation whose correction to scaling behaviour is 
given by ~ L a ~ d . For Tp = T c , the role of the two terms in the right hand 
side of (3.23) is reversed. The thermal contribution provides the dominant 
scaling term ~ L a while the first term yields correction to scaling ~ L a ~ d . 
Taking A = and setting the coupling constant g at the fixed point value 14 ' 15 
the first term can be made to vanish (t — > 0) obtaining 

C(k,t) ~ L a T c F € (x) (3.24) 

with 

F e (x)= f dy(l-y)-'l a e- wxay . (3.25) 
io 
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Notice that lim e ^o F € (x) = Fq(x). Thus, at the upper critical dimensionality 
2a the non trivial fixed point merges with the trivial fixed point on the 
temperature axis. 

In summary, with NCOP we have found the following asymptotic scaling 
properties 
1) [7> = 0,Aii] 

C(k,t) ~L\t)F{x) (3.26) 



2) [T F = T C > 0,/n] 

3) [T F = 0,/i 2 ] 



C(k,t) ~ L°{t)T c F Q {x) (3.27) 



L e (t)F(x) for d > d c 
C(M)~< for d = d c = 9 + a (3.28) 



log £,(*)' 

L d - a (t)F{x) for d < d c 



4) [T F = T c >0,» 3 ] 



% \ L°(t)T c F (x) for d > 2a 

j L°(t)T c F £ (x) ford<2a (3 ' 29) 

5) [T F <T C ,H 

C(k,t) ~ L d (t)F(x) (3.30) 
where L(t) ~ t 1/,T , a; = kL(t) and 

F(x) = e-^// (3.31) 

F (x) = (l-e- wx °)/wx a (3.32) 

FJx) = C dy(\-y)- e l°e- wx ° v (3.33) 
Jo 

with e — 2a — d. 
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3.2 - COP 



Let us now go back to Eq(3.1) with p ^ 0. The important novelty is that 
now by dimensional analysis we can form the two lengths L(t) = (2Tt)~ 
and X(t) = (2T \ Q \) 1 ^ p . In order to establish the scaling behaviour we must 
know which of the two is the dominant one. 

T F = 

Considering first, as before, quenches to zero temperature the structure 
factor is given by 

where a = sign(Q). To begin with let us rescale lengths with respect to L{t) 

C{k, t) = Ae- af3xP L e F > (x) (3.35) 

where 

P=(jY (3-36) 
P >( x ) = —JT- (3-37) 

Jb 

with x = kL(t) and the reason for the notation F>(x) will be clear below. 
Again, we must solve for Q(t) using the analogue of Eq.(3.4) 



^ = r + gAL e - d K d dxx^e-^^^. (3.38) 
at Jo 



l i 1 = {r = 0,g = 0) 
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Since in this case Q(t) = implying \(t) = 0, clearly L(t) is the dominant 
length and z = p + a. From (3.35) we have, as in the NCOP case, that 
standard scaling is exactly obeyed over the entire time history with 

C(k,t) = A^F^x). (3.39) 

When the nonlinearity is present it is no longer possible, contrary to 
NCOP case, to solve directly for Q(t). Hence, we now make statements about 
the solution of (3.38) for large time by consistency checks on the assumption 
that either one of the two lengths L(t) and \(t) prevails over the other. 



lii = (r = 0, g > 0) 



From (2.5) and (2.8) in this case a is positive. Let us first suppose that L 
prevails over A, i.e. f3(t) — > as t — > oo . For large time Eq.(3.38) can be 
replaced by 

^ = gAL e ~ d K d r dxx^e-^ (3.40) 
dt Jo 

and integrating we find 

Q(t) ~£^- d +mnst (3 - 4i) 



where d c = 6 + a . Inserting this result into (3.36) follows 
/3(f) ~ < 



' L- p for d > d c + p 

L-p log L ioid = d c +p (3.42) 
L d c -d for d<d c +p 



which is consistent with the assumption if d > d c . Thus, for d > d c from 
(3.35) we find 

C{k,t) ^ L e [l- p{t)x p ]F > (x). (3.43) 
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This result together with (3.42) shows that there exists yet another critical 
dimensionality d c = d c + p affecting the behaviour of corrections to scaling. 

Assuming next that (3{t) takes a constant value c, we find that this is 
consistent only for d = d c and from (3.35) follows 

C(k,t) ~ L 9 F{x) (3.44) 

with 

F(x) = _L e -[- p+CT +- p ] (3.45) 

showing that there is scaling with the same exponents as at ^ , but with a 
modified scaling function. 

Finally, assuming (3(t) — > oo , from (3.38) follows 

(3{t) ~ LT^+p) (dc ' d) (3.46) 

which for consistency requires d < d c . Therefore, now we must scale with 
respect to A obtaining 

C(k,t) ~ X e [l - wp-'i^x'^F^x') (3.47) 

where 

F<(*') = ^ (3.48) 
and x' = k\(t), X(t) ~ t l ' z with 

z = d + p-9. (3.49) 

This result is interesting because a pattern qualitatively different from the 
corresponding case with NCOP is obtained. In the latter case when the 
nonlinearity becomes relevant below d c the exponent a changes from the 
value 6 to the dimensionality dependent value (d — a) , while the exponent z 



20 



and the scaling function remain the same as for the quench to the trivial fixed 
point /ii. With COP instead we find that a keeps always the trivial value 0, 
while there is a change in the scaling function and the growth exponent picks 
up the value (3.49) dependent on the space dimensionality of the system. 



Ms = (r < 0, g > 0) 



For the quench in the phase ordering region Eq.(3.38) can be rewritten in 
the form 

f = r + #»w> < 3 - 50 > 

where 

POO 

I(f3) = / dxx^^e-MW (3.51) 
J o 

f(x) = ax p + w/T V +(T . (3.52) 

Making the assumptions to be verified a posteriori that (3 is asymptotically 
divergent and that a is negative we can make a saddle point evaluation of 
(3.51) 

7T 



where 



and 



u 



wa(p + a) 



-U -v 

e p u 



7 = [2(d 



w(p + a) 



(p + <j)]/2(p + <j). 



(3.53) 



(3.54) 



(3.55) 



For large time we can set to zero the left hand side of (3.50) obtaining 1(0) = 
—rL d ~ e /gK d A . Next, using (3.53) and taking the logarithm we find the 
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asymptotic relation 



p(d -6) py 

u— logL Xogu. (3.56) 

wa wa 



Inserting the leading contribution into (3.54) we obtain 

(3(t) ~ (logL)(^T (3.57) 

consistently with the assumption. Thus, in this case X(t) and L(t) diverge in 
the same way up to a logarithmic factor. As fundamental length we choose 8 
the inverse of the wave vector k m (t) where C(k,t) reaches its maximum 
value 

k m (t) = L -1 upfe. (3.58) 

Inserting this result in (3.35) and introducing the variable x = k/k m we 
obtain 

where 

(p(x) = V -^* v - x p+a - (3.60) 

Next iterating once (3.56) and inserting the result into (3.59) we obtain the 
asymptotic expression 8 ' 10 

~kt d {k m L)^ 



C ^-~ fejyT ( 3 - 61 ) 

which, up to a logarithmic factor in the amplitude, is in the multiscaling 
form (1.5) 

C(k,t)~k-^± (3.62) 



with 



a (x) = — — —(p(x) + 9. (3.63) 
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This exponent ao(x) is plotted in Fig. 6 for 9 = 0, a = 2 showing that the 
transition from multiscaling to standard scaling is a smooth one as p is varied 
continously from p = 2 to p = 0. 

Notice that from (3.58) and A;" 1 ~ t 1 ^, the exponent z develops a time 
dependence which asymptotically is given by 

(3.64) 
T F > 

fii = (r = 0,g = 0) 

Again, with Q(t) = and 9 = it is straightforward to derive from (3.1) the 
standard scaling result 

C(k, t) = Ae- wxP+a + L a T F F (x) (3.65) 

where 

F (x) = (1 - e- wxP+a )/wx a (3.66) 
and the same considerations made about (3.19) apply here. 

Ms = (r < 0,g > 0) 

For quenches to /i 3 = (r < 0, g > 0) now we cannot solve explicitely for 
Q(t). Then we proceed differently by rewriting (3.1) as 

C(k, t) = Ae u( * Mi) + 2FT F k p C ^'e^M^-^'M*') (3.57) 

Jo 



z ~ (p + a) 



1 + 



log log t 
logt 
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where (p(x) has been defined in (3.60) and u(t) is related to Q(t) by (3.54). 
In the following we will drop the bar over x. If by analogy with what we 
have found at T F — we make the ansatz 

e «(t) = cL p^ ( 368 ) 

with p > 0, the quantities to be determined are the constant c and the expo- 
nent p. Inserting into (3.67), and assuming k^{t) ~ L(t) up to logarithmic 
factors we find 



C(k,t) = A(cL p ) 



ip(x) 



>) e -p V (x>)logL 



l+T F L a — ( X dx'x'^-Ux/x'Y^ 
wx a Jo 

(3.69) 

For large L the integral in the right hand side can be evaluated by steepest 
descent. Keeping in mind that p > and that <p(x) behaves as in Fig. 6, we 
must distinguish between x < x* and x > x* where x* = ( 2 ^ £ ) 1 ^ cr is the non 
trivial zero of (p(x). For x < x* the exponential reaches the maximum value 
at x' — 0, while for x > x* the maximum is at x' — x yielding 

C(k, t) ~ A(cL p yW + ^^{cL P Y {x) (3.70) 

for x < x* and 

C(k, t) ~ A(cL P Y {x) + (3.71) 

fro x > x*. 

We now determine c and p by using the self-consistency condition (2.14). 
From (3.54), where f3 is defined in (3.36), and (3.68) follows R(t) = Q(t) ~ 

a 

(plogL)p+ CT L~ a . On the other hand R(t) can be computed from (2.5) with 
S(t) obtained by integration of (3.70) and (3.71) over k. Thus, for T F = 
dropping logarithmic factors one finds 

rAL 

L-° = r + gAK d L~ d / dxx d ~ l e^°^ cLP \ (3.72) 
Jo 
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For large L the integral picks up the dominant contribution at x — 1 giving 

L- a = r + gAK d c a/p L {pa - pd)/p . (3.73) 

The right hand side can vanish only if p = pd/a, in agreement with (3.63), 
and c CT /f ~ Ml 

For Tp 7^ 0, the first terms in the right hand side of (3.70) and (3.71) are 
asymptotically negligible with respect to the second ones yielding in place of 
(3.72) 

rAL 



L-° = r + gT F L°- d K d 



/ dxx d - a -\cL P Y {x) + / dxx d ~ 

JO J x* 



•o— 1 



(3.74) 



Furthermore, keeping into account that for large L the dominant contribution 
to the first integral is obtained at x — 1 and that the contribution at x* is 
negligeble in the second one finally we obtain 

L~° = r(^— ^) + gT F K d c°l p L p ^°- d . (3.75) 

J- c 

For Tp < T c this implies p = p(d - a) /a and c a / p ~ M 2 (^^). Inserting 
into (3.70) and (3.71) one finds 

C '(*'*)~T>— (3.76) 

with 

/ x f ~(d — a)(p(x) + a for x < x* , n __ x 

a(x) = < a r * (3.77) 

w [ a for x > x* v ' 

showing that the structure factor obeys multiscaling for x < x* and standard 

scaling for x > x*. Notice that the two behaviours match at x* since (f(x*) = 

0. 

If the quench is made on the critical surface (T F = T c ), then Eq. (3.75) 
gives p = £(d — 2a) which is negative below the upper critical dimensionality 
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contradicting the initial assumption p > 0. Hence we take p = which 
implies standard scaling with u(t) constant. From (3.68) and (3.36) we have 
Q(t) = bL p {t). Inserting in (3.1) and taking A = for simplicity, we obtain 

C(k, t) = L a T c F e (x) (3.78) 

with 

FJx) = — [ X dx ' x 'p+°-i e -l^ + °-*' p+ n+b^-z'n]. (3.79) 

K ' wx u Jo v ' 

From the self-consistency condition (2.14) now we find 



2Yp 



• ( r°° j too ^ 1 "i 

-gK d I?°- d T c / dxx d -'[F e {x) - — ] - / [F e (x) - — ] 

UO «JX ff JAL HJX J 

(3.80) 

where we have used the definition (2.16) of T c . From (3.79) for large x we 
have F e ~ ^(1 - cfcr- CT ) where c = ^ J °° d^e'^ . Inserting into (3.80) 
and defining /(&) = / °° tte^-^F^a:) - -M we have 



2r P 



gK d L^ d f(b) (3.81) 



2rp(rf - 2(t) 

from which follows b = for d > 2a, while for d < 2a the value of 6 is given 
by the condition f(b) = 0. Solving this equation in the e-expansion 15 one 
finds b ~ e as for NCOP, which implies the same structure of fixed points. 
Namely, for e — > the non trivial fixed point still merges with the trivial one 
and lim^o F e (x) = F (x). 

The summary of the asymptotic properties with COP is given hereafter 

1) [7> = 0,^i] 

C(k,t) ~ L e (t)F(x) (3.82) 

2) [T F =T C)/Ul ] 

C(M) ~ L a (t)T c F (x) (3.83) 
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3) [T F = 0,fi 2 ] 



L e (t)F>(:r) for d >d c 
C{k, t) ~ { L e (t)F(x) for d = d c = + a 
\ 9 (t)F < (x') iord<d c 



where x = kL(t) and x' = k\(t) 
4) [T F = T c >0,H 



c(M) 

5) [0<T F <T c?y u 3 ] 



L a (t)T c F Q {x) for d > 2a 
L°(t)T c F £ (x) iord<2a 



C(k,t)=T F - ^ 

v ' wx a 



6)[T F = 0,fi 3 ] 

C(k,t) ~ L Q0 ^(t)^ 



with 



X(t) ~ iV^-f) 



x' e 



F (x) = -±-(1 - e~™ P+ff ) 



F e (x) = — ( X dx'x'^^e-^^-^^ 



a (x) = — — — +0. 



2 + (d - 2)v?(x) for rr < a;* 
2 for rr > a;* 
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(p(x) = ^—^x p - x p+a . (3.97) 
V 

4 - Renormalization group 

We now discuss the scaling properties of the model within the RG approach 
to the problem 16 . Let us recall that in static critical phenomena the Wil- 
son RG equations are obtained performing the following operations on the 
equilibrium probability distribution P eq [(j); Tp, //] 

i) elimination of hard modes (f>(k) with A/ 1 < k < A where / > 1 

ii) rescaling of wave vectors and order parameter 

k' = lk 

= l-*$(h) 

iii) requirement of form invariance of P eq [4>', T F , fj] . 

These operations generate recursion relations for the parameters (Tp,/x) 
which allow to describe scaling in terms of the geometry of the fixed points 
and their domains of attraction. 

In quench processes one deals with a time dependent probability distribu- 
tion P[0;t,Tp, jj] . Therefore, RG transformations performed on this object 
are expected to give recursion relations for (t, Tp, /i) with a fixed point 
structure which accounts for the variety of scaling behaviours obtained in 
section 3. In order to implement the procedure outlined above, we should 
construct P[<j>; t, Tp, fj] and then carry out renormalization. However, since 
the stochastic process is gaussian and all the equal time information is in 

— * 

C(k,t) , we can work directly with the equation of motion (2.13). First we 
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separate soft and hard modes 

C(k,t) = C s (k,t) + C h (k,t) (4.1) 

with 

C(kt)-\ C ^ t ) forO<A;<A// 

MM)-| for A// < k < A (4 ' 2) 

,p , ( for < k < A/l 

The equation of motion for either component is given by 

dCs ^^ = -2i>P+ CT + k p R(t)]C S;h (k, t) + 2Tk p Tp (4.4) 

with 

R(t)=r + gS s (t)+gS h (t) (4.5) 

S s , h (t)= J -0- d C s , h (k,t). (4.6) 
Then we proceed to eliminate hard modes. Integrating the equation of motion 

— * 

for C h (k, t) we find 

C h (k,t) = C h (k,0)e- 2 ^ wkP+at+kP ^+2Tk p T F C dt > e -^[^(t-t') + ^m)-Q(t'))] 

Jo 

(4.7) 

— * 

where Q(t) = Q s (t) + Qh(t) ■ After integrating over k we should solve 
for Sh(t) and insert the result into the equation for C s (k,t) . However, 
considering that we must eliminate modes with k > and that in the 

scaling regime these have already equilibrated, with a good approximation 
we can set 

c, ( M)^ A cc) = — (4 ' 8) 

— * 

and we are left with Eq.(4.4) for C s (k,t) with 

R{t)=r + gS a (t)+gS h {oo). (4.9) 
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Next we carry out the rescalings for k < A/7 

C'(k', 0) = l d ~ 2y ^C(k, 0) 

k' = Ik 
t' = r z t 

C'(k',t') = l d - 2y(x) C(k,t) 



(4.10) 

(4.11) 
(4.12) 
(4.13) 



where we have allowed for an order parameter scaling index dependent on 
the invariant quantity x = kt x l z . From (4.4) and (4.13) we obtain the 
transformed equation of motion 



dC'(k',t') 



-2F 



dt' 

+2Tk' p l z - p+d - 2y{x) T F 



l z - a - v wk' v+a + l z - p k' p R(l z t') + r" 1 ^^ log/ 



C'(k',t') 
(4.14) 



which can be rewritten in the same form as the original equation of motion 
(2.13) 



f}C'(P t') 

{ -^- = -2T[w f k' p+(7 + k' p R'(t')]C\k', t 1 ) + 2Yk' p T' F (4.15) 



dt' 



defining 



w ' = l z ~°- p w 

r^y loe-/ 
R\t') = l z - p R(l z t') + d * 7 

J- Ai v 

rpf^ _ ^z-p+d-2y(x)rp^ 

In order to preserve the self-consistent structure we must require 

R'(t') =r' + g'S'(t') 
d d k' 



(4.16) 

(4.17) 
(4.18) 



S'(t>) 



(2tt)< 



C'(k',t'). 
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Since the left hand side of (4.17) depends only on t' , the x-dependence on 
the right hand side must disappear implying 



This form of y(x) diverges at x = and we must necessarily have c = 0. 

So far we have managed to map the process described by (2.13) in the 
new process governed by (4.15) with new parameters A', w', //, T' F . We are 
going to be interested in those processes (fixed points) whose parameters 
A*,w*, n*,Tp do not change under renormalization. For these values of the 
parameters, since the form of the original equation of motion (4.15) and the 
self-consistent structure have been preserved we have C'(k',t f ) = C(k, f ,t'). 
The existence of such fixed points is compatible only with the standard scal- 
ing choice y(x) = y, which seems to exclude the multiscaling solution found 
in the previous section. We shall comment on this later on. 

We proceed to extract recursion relations. From (1.4) and (4.10) follows 



y(x) = clog a; + y(0). 



(4.19) 



C'(k', 0) = i d+e - 2y A/k 



(4.20) 



implying 



A' = i d + e ~ 2 v/^ 



(4.21) 



Next, rewrite (4.17) as 



r' + g'S'it') = l z " p [r + gS s {l z t') + gS h (oo)} 



(4.22) 



where for quenches to final states with R(oo) = one has 



S h (oo) = T F B(0)(l-r- d ) 



(4.23) 



and -B(O) is defined in (2.18). Since from (4.13) we have 



S'(t') = l 2id - y) S s {t) 



(4.24) 
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inserting (4.23) and (4.24) in (4.22) we find r' = l z ~ p [r + gT F B(0)(l - 1^)} 
and g' = l z -p+ 2 (y- d ) g . Introducing the scaling field 

r = r + gT F B(0) = r(T c - T F )/T C (4.25) 

the whole set of recursion relations is given by 

A' = l d+t) - 2 yA (4.26) 

w ' = l z ~°~ p w (4.27) 

t' = l z ~ p r (4.28) 

g' = l z ~P+ 2 (y- d ) g (4.29) 

T' F = l z - p+d ~ 2y T F . (4.30) 

We emphasize that the use of the same scaling index y for t > as well as for 
t — leads to fixed points as processes for which scaling holds over the entire 
time history. For processes where scaling invariance is only asymptotic is not 
necessary that there exists a non trivial fixed point solution of (4.26-4.30). 

Fixed points 

The next step is to look for fixed points (A*, w*, r*, g*, T F ) of the recursion 
relations and to extract from them the exponents z and a. These exponents 
are determined by imposing that two of the parameters have a finite fixed 
point value. In the remaining recursion relations we consider the trivial fixed 
point solution and the corresponding domain of attraction. In order to show 
how this works in practice let us begin by imposing that A* and w* be finite. 
From this, using (4.26) and (4.27) follows 

z=p + a (4.31) 
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and 

a = 2y - d = 6. 
Inserting these values into (4.28),(4.29)and (4.30) we obtain 



(4.32) 



t' = Ft (4.33) 

g' = l dc - d g (4.34) 
T' F = r~ e T F (4.35) 
with d c — 6 + a. The trivial solution is r* = g* — Tp — which coincides 



with (Tp = 0,/ii). The corresponding domain of attraction, considering that 
for quenches from high temperature to the critical surface 6 < a, is given 
by the g-axis, i.e. (Tp = 0,/i2) for d > d c . Otherwise, ford < d c , this fixed 
point is completely unstable, and (4. 31), (4. 32) do not apply for quenches to 
(Tp = 0,/i 2 )- This coincides with what we have found in section 3. 

Next, let us impose that A* and g* be finite. Then we find z = d + p — 9, 
a = 9 and the remaining recursion relations 

w ' = l d ~ d c w (4.36) 

t' = l d - e r (4.37) 

Tp = l d - 2d T F (4.38) 

with the solution w* — t* — Tp — 0. w flows to zero for d < d c . This fixed 
point for COP corresponds to the quench to (Tp = 0,/i 2 )- For NCOP this 
corresponds to the case of independent particles and it is possible to show 
that in this case the amplitude of the structure factor vanishes. Thus in order 
to treat quenches to (Tp = 0, /x 2 ) for NCOP with d < d c we must require 
that w* and g* be simultaneously finite. This reproduces immediately the 



33 



results z — a and a = d — a of (3.13). The ensuing recursion relations for 
the other parameters 

A' = &- d A (4.39) 

t' = Ft (4.40) 

T' F = l 2a - d T F . (4.41) 

yield a trivial solution which is unstable under all perturbations. The mean- 
ing of A flowing to infinity can be understood from the result of section 3 
where the crossover time t* vanishes for A — > oo. In this limit we have 
a fixed point in the sense specified above that the same scaling behaviour 
applies over the whole history of the process. 

So far we have dealt with fixed points with T F — and r* = 0, namely 
with quenches to zero temperature critical points. In order to analyse quenches 
on the critical surface at finite temperature, we must require T F and w* fi- 
nite as it is usually done in static critical phenomena. From these conditions 
follows z = p + a, a = a and 

A' = l e - a A (4.42) 

t' = Ft (4.43) 
9 = l 2a - d 9- (4.44) 

Thus, for quenches to finite temperature on the critical surface A is irrel- 
evant and the attractive fixed point goes from trivial to non trivial as the 
dimensionality goes from above to below the upper critical dimensionality 
2a. 
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Finally let us come to the discussion of quenches inside the phase ordering 
region. For this we require that w* and the fixed point ratio 



.9, 



M, 



0V T r 



(4.45) 



be finite obtaining z = p + a and a — d. The other recursion relations 

A' = l e - d A (4.46) 

T' F = r~ d T F (4.47) 

show that temperature perturbations in the ordering region are irrelevant 
and that A flows to zero. Now, with A* = Tp = the structure factor 
C(k,t) vanishes identically, namely the scaling form (3.18) is obeyed with 
F(x) = 0. This means that a non trivial scaling solution for a quench in 
the phase ordering region with a = d is necessarily asymptotic and cannot 
be made to hold over the entire history of the process by any choice of the 
parameters. 

This completes the analysis of the fixed point structure of the phase 
diagram and the derivation of exponents. The RG treatment of the problem 
presented above reproduces the whole structure found in section 3 except 
for multiscaling in the quench below T c with p ^ 0. The point is that the 
RG procedure we have followed above yields the exponents z and a within a 
standard scaling framework, but gives no information on the scaling function. 
If one goes further by performing the scaling ansatz in the equation of motion 
an equation for F(x) is obtained and it turns out that the scaling function 
vanishes if T F < T c and p ^ 0. In order to recover multiscaling trough the 
RG approach the set of transformations must be properly generalized. We 
do this only for the (T F = 0,/i 3 ) case 17 . 
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As we have seen in section 3, when there is multiscaling z is weakly time 
dependent. Let us then generalize the set of transformations (4.10-4.13) by 
allowing z to depend on t with the constraint 



lim z(t) = z^ = p + a. 



t^oo 



(4.48) 



With these modifications and T F = in place of (4.14) we find 

rid' (P t'\ 

{ - ' ' = -2T[wk' p+ ° + k' p FR(t)}C'(k',t') 
of 

-2 \y(z - z 00 )[wk ,p+ ' 7 + k' p l a R(t)\ + j C'(k',t') log/. (4.49) 

Imposing the requirement of form invariance we find that (4.49) is of the 
form (4.15) if in place of (4.17) the following conditions are satisfied 



dy 

dt 1 



R\t') = FR{t) 

-r[ z (t') - 2co]K' p+a + k' p R\t% 



This latter equation holds also for unprimed variables and using 



dx x 
dt zt 

we obtain 

dy Tz[z - z^t 1 '^^'* 



t dz , 
l--—\ogt 

z dt 



dx 



[1- If l°g*l 



wx 



.P+cr-1 , p-l.a/z 



+ x p - 1 t a/z R(t)) . 



(4.50) 
(4.51) 

(4.52) 



(4.53) 



In order to get rid of the time dependence on the right hand side we must 
have 

t dz 

- Tz(z - Zoo ) = S z - Z ^l z {\ - -— log*) (4.54) 

z dt 

t°l z R(t) = b (4.55) 
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where b and c are constants. Integrating (4.53) we find 

r p+a £ 



y(x) 



WX' J 

h -x p 

p + a p 



+ 1/(0). 



(4.56) 



From (4.55) and the definition (2.5) of R(t) follows that the sign of b is 
determined by the parameters [i = (r, g) . In particular b is a negative 
quantity at fi^. Imposing y max — d , as appropriate for quenches to 113 , 
taking the position of the maximum at x — 1 and using y(0) = (d + 0)/2 
from (4.56) we find b = —w , c = {6 — d)p(p + a) /2a and 



2y(x) — d = a (x) 



(4.57) 



where a (x) coincides with (3.63). Hence, for quenches to /i 3 , by allowing 
for a time dependence in the growth exponent z , we have recovered via RG 
the multiscaling behaviour of the exact analytical solution. 

The explicit time dependence of z(t) is obtained extracting the asymptotic 
behaviour from (4.54) 

log log t 



z(t) ~ z Q 



1 + 



log* 



(4.58) 



which is consistent with the assumption (4.48) and reproduces (3.64). Finally, 
inserting this result into (4.50) we find 



R(t) 



log A p+a 



(4.59) 



in agreement with (3.57). 
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5 - Conclusions 



In this paper we have investigated in detail the solution of the large-N model 
for growth kinetics with the aim of giving a comprehensive view of the in- 
fluence on the scaling properties of the various elements which enter in the 
specification of the problem. These are the presence or the absence of a con- 
straint on the order parameter (COP or NCOP), the initial condition, the 
structure of the phase diagram of final equilibrium states, the range of the 
interaction and the dimensionality of space. 

What the model shows, apart for quenches to the trivial fixed point, is 
that scaling properties are quite different with and without conservation law. 
This is due to the existence of only one divergent length L(t) for NCOP and 
of two divergent lengths L(t) and X(t) for COP. It is the interplay between 
these two lengths which leads to phenomena not observed with NCOP such 
as i) a change in the growth law when crossing the critical dimensionality 
for the quenches to (Tp = 0,7x2) and ii) multiscaling for quenches inside the 
phase ordering region. About multiscaling, the availability of the rich variety 
of cases illustrated in the paper should allow to speculate about its origin. 
Thus we have found that for NCOP in no circumstance there is multiscaling 
and the same holds true for COP, except for the quenches below T c . What, 
then, is the peculiarity of these processes. One possible interpretation is that 
in these processes the system orders and likes to do so by condensing, i.e. 
by growing a peak which scales like L d , at k — 0. This is fine with NCOP, 
but with COP there is a conflict with the conservation law which prevents 
anything to happen at k — 0. In this case the peak is formed at some k m 7^ 0. 
The compromise realized in the large-N model is multiscaling whereby the 
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behaviour L a °^ of the structure factor interpolates smoothly between the 
behaviours LP at k = and L d at k m . This picture fits nicely with the absence 
of multiscaling in any process with NCOP and in all processes with COP on 
the critical surface. In fact, in the latter case there is nothing to condense at 

— * — * 

k = and in the former there are no constraints at k — 0. However, with this 
mechanism for multiscaling there should be nothing special about N = oo 
and multiscaling should be found also for iV < oo. Numerical simulations 18 
so far have reported no evidence for multiscaling for N = 1 and iV = 2 in 
two and three dimensions. This could mean that multiscaling disappears for 
N < d, when localized topological defects appear in the system. However an 
argument against the existance of multiscaling for any finite N is the result of 
Bray and Humayun 19 . By analysing an equation of motion for the structure 
factor which includes the first order correction in l/N, they have reached the 
conclusion that multiscaling does not survive for iV < oo since the correction 
term sustains a standard scaling solution. 

What we can say is that both standard scaling and multiscaling imply 
scale invariance of the structure factor due to the presence of a divergent 
length. So far we have not found a criterion to predict a priori which one 
should hold. Only a direct calculation, either analytical or numerical can dis- 
criminate between the two. It must be emphasised that these concepts apply 
also to other models. For example in DLA (diffusion limited aggregation) 
model in two dimensions it has been found numerically 20 that multiscaling 
holds. In any case, even if there is not a general criterion, multiscaling should 
be more likely to occur in situations where the width of the interface becomes 
very large. 

Let us then comment on those features of the N = oo solution which we 
believe to be of general validity. The crossover structure which emerges as 
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the parameters of the quench are moved over the manifold of final equilib- 
rium states is a generic feature which is expected to hold beyond the large-N 
model. The main point of our analysis is that it is quite possible, before 
the true asymptotic behaviour is reached, to detect a preasymptotic scal- 
ing behaviour due to a less stable fixed point lying in the neighbourhood of 
the final equilibrium state. It is natural to pose the question of the observ- 
ability of these effects. Here we suggest (Appendix I) that the duration of 
this preasymptotic behaviour can be magnified and observed in off-critical 
quenches 21 . 

Furthermore, the crossover picture we have illustrated suggests the pos- 
sibility of observing a crossover in the growth law (1.2) in the symmetrical 
quench of a system with scalar (N = 1) COP. In that case asymptotically 
L(t) grows according to (1.2) with z — 3. On the other hand in the triv- 
ial theory (r = 0,g — 0) one has z = 4 for COP, irrespective of the order 
parameter being a scalar or a vector. Thus, for a quench to (T F = 0,/i 3 ) 
sufficiently close to (T F = 0,/ii) with TV = 1 it should be possible to observe 
the influence of the trivial fixed point at early time producing a crossover in 
(1.2) from z = 4 to z = 3. 

6- Appendix I 

Let us consider a process where symmetry breaking along one direction is 
allowed, e.g. 

<<P a (x,t)>=N 1 / 2 M(t)5 a>1 . (6.1) 
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This may be due to non symmetrical initial conditions, or to the presence 
of an external field or to both of these circumstances 22 . Introducing the 
external field along the 1-direction from (1.3) dropping the long range term 
we obtain the equation of motion for the order parameter 



d<p a (x,t) 
dt 

+T] a (x, t) 



-r(?v) p 



-V 2 + r + |p ]T N4>1 (x , t) ] Q (f, t) - K 



(6.2) 



where h a = N l l 2 h5 a ^i . 

Defining the fluctuation field ip(x,t) by 

<f>x{x,t) = N^ 2 M{t)+^{x,t) 

and inserting into (6.2) we obtain the pair of equations 

d{N l l 2 M + V>) 



(6.3) 



dt 



-T(iV) p -VV + r(iV 1/2 M + ^) + 



jjiN^M + V) E * + ^(iV 1/2 M 3 + 3M 2 ^ + ^M^ 2 + l^ 3 ) - 



+ T) 1 (x,t) 



(6.4) 



dt 

+— i/> 



-r(iv) p 



-V 2 0« + rc 



7^1 



(6.5) 



Taking the large-N limit we replace ^^(a 7 , i) by < t) >— S±(t) 

and collecting terms of the same order of magnitude, from (6.4) and (6.5) we 
obtain the set of equations 

= _r(*V) p ([r + <?M 2 (t)]M(t) + gS ± (t)M(t) - h) (6.6) 
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^M. = -r(<V)* [-V 2 + r + 3gM\t) + i>(x, t) + ^(f , t) (6.7) 

dM d f t] = -r(<V)' [-V 2 + r + ^M 2 (t) + <?S ± (i)] M2>t)+Vp(2,t). (6.8) 

Notice that the behaviour of ^(x, t) is immediately obtained once the pair of 
coupled equations (6.6) and (6.8) have been solved. Furthermore, for COP 
Eq.(6.6) is trivial since M(t) does not change in time and keeps the initial 
value M(0) . In this case the equation (6.8) for the transverse components, 
after Fourier transforming can be rewritten as 

d(f) (k,t) 



dt 



= -Tk p 



-V 2 + f + gS ± (t) (/> (x,t)+r] (x,t) (6.9) 



where 



f = r + gM(0) 2 . (6.10) 
Thus f can be modulated by varying M(0). In particular a quench to /i2 



can be realized as an off-critical quench to Tp 
coexistence region M 2 (0) = —r/g. 



and at the edge of the 



7- Appendix II 

In the limit t — > oo the left hand side of Eq.(2.13) vanishes and we have 

= -2V[wk p+a + k p R]C{k, oo) + 2Tk p T F (7.1) 

where R stands for R(oo). 

Let us first consider p = and a system in a finite volume V . Assuming 
R > from (7.1) follows 
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where R must satisfy the self-consistency condition 



fl = r + J-V 7 Tf (7.3) 

k 

which always admits a solution with R > . 

In the infinite volume limit Eq.(7.3) can be rewritten as 

R = r + gT F B(R)+g^ (7.4) 

where, allowing for the possibility that the solution R — > , the zero wave 
vector term in the summation has been separated out and 

B (R) = [ d \ - 1 - (7.5) 

is a non-negative monotonously decreasing function of R with a maxi- 
mum at B(0) = K d A d ~ a /[w(d — a)] . A is a momentum cutoff and K d = 
(2 d - 1 ir d / 2 T(d/2))~ 1 . From Eq.(7.4) follows that when r < there is a critical 
temperature 

T c = ^— (7.6) 

gB(0) 1 J 

such that for Tp > T c there exists a solution with R > and therefore the 
last term in the right hand side can be neglected, but for T F < T c Eq.(7.4) 
can only have the solution R = provided that R vanishes in the infinite 
volume limit as it! ~ 1/V . Hence, defining the constant 

M* = H (77) 

and inserting into (7.4) we find 

M 2 = Ml{T c -T F )/T c (7.8) 
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with Mq = —r/g. Notice that for Tp = T c Eq.(7.4) admits the solution 
R = provided R ~ V~ x with < rr < 1 . In conclusion the structure 
factor is given by 

r<h \- \ T F^ wk<J + R ) with R > for t f > T c ( 7 Q\ 

[ ,OC) -\T F /wk° + (27r) d M 2 5(k) for Tp <T C . ^ 

For p ^ Eq.(7.2) applies only for k ^ 0. Due to the conservation law 
C(k — 0, oo) is determined by the initial condition. If we consider an initial 

— * — * 

state without symmetry breaking we have <f>(k = 0) = and C(k — 0,t — 
0) = 0. Then Eq.(7.3) must be replaced by 

R = r + ^ 7 T / , Tf - . (7.10) 
V wk a + R y ' 

Solving the above equation for R we find that there exists a temperature 

f(V) = -- — (7.11) 

such that R > for T F > T(V), while R < for T F < T(V) with \R\ < 
k min where k min ~ V~ x l d is the minimum value of the wave vector in the 
summation. When the infinite volume limit is taken from (7.11) follows 
f(V) -> T c and the anlogue of (7.4) is 

R = r + 9 T F B {R) + ^-^-- R . (7,2, 

For Tp > T c there is a solution R > and the last term vanishes, while for 
Tp < T c Eq.(7.12) admits the solution R = provided (£;^ n + i?) ~ l/V. 
Writing 

M 2 = — , Tf — (7.13) 

V(iuA& in + i2) 

in place of (7.7), in the end we recover the results (7.8) and (7.9). 
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Figure captions 

Fig.l - Renormalization group flow on the temperature axis. 

Fig. 2 - Manifold of final equilibrium states with the critical surface separating 

disordered states (above) from ordered states (below). 

Fig. 3 - Behaviour of S(t) in a quench to /i 2 for NCOP with a = 2, A = 10 
and d > d c : (d — 3,9 — 0). The straight dashed lines have slope —3/2. 
Fig.4 - Behaviour of S(t) in a quench to fi 2 for NCOP with a — 2, A = 10 
and d < d c : (d = 3, 9 = 2). The top dashed line has slope —0.5, the others 
below have slope —1.0. 

Fig. 5 - Behaviour of S(t) in a quench to /i 3 for NCOP with a — 2, A = 
10, d = 3,9 = at fixed g = 1 as r approaches the /x 2 -axis. The straight 
dashed line has slope —1.45. 

Fig. 6 - Plot of a (x) for different values of p. Multiscaling goes over smoothly 
to standard scaling as p — > 0. 
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